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1. Introduction 



Two dimensional quantum gravity has been solved non-perturbatively by taking 
the double-scaling limit of the matrix model - the dual model for the 2d discretized 
gravity. Here one has a parameter which determines the way this double limit 
is approached. The potential is fine-tuned as the size of matrix gets large, and the 
contributions from difi^erent 2d topologies may be obtained as a power series in 
with the power of t<~'^ being the Euler number of the Riemann surface. Surprisingly, 
the scaling property (dependence on the cosmological constant) of the expansion 
agrees order by order with the genus expansion in the 2d quantum gravity coupled 
with conformal matters.* This is the reason why k is called the string coupling 
constant. Furthermore, one can determine the scaling limit non-perturbatively by 
solving a characteristic non-linear ordinary differential equation (string equation) 
satisfied by the string susceptibility. 

Even though one can find infinitely many critical points in the variety of 
models and their potentials, which all have double scaling limits, these continuum 
limits can be classified into hierarchical universality classes due to the integrability 
emerging at these critical points. Actually integrability is there even at the discrete 
level. If the size of matrix is finite, say N ^ the hermitian one matrix model can be 
embedded into the GL{N) Toda system, with variation of the coupling constants 
of the matrix model forming the Toda flows [3] , so that the continuum integrability 
may be understood by the scaling from the discrete integrability. 

1.1 GL{N) Toda Hierarchy 

We will explain the Toda flow briefly by following refs. [1,3] and Neuberger 
[4]. After integrating out the degrees of freedom for the similarity transformations, 
one may evaluate the matrix integral 

Z=^ f d^Vexp[-/?trF(0)] = (vac|vac), (1.1) 

-'^ • J NxN hermitian 

in terms of the A'^-fermion vacuum 

Po(Ai) ••• Po(Aiv) 



vac ) = 



/N! 



Piv-i(Ai) ••• Pn-i{Xn) 



(1.2) 



where the (n -|- l)*'^ state Pn(A) is an orthogonal monic polynomial of degree n 
deflned by 

/oo 
dX e-^^(^) Pn{X) Pm{X) = hn 5n,m. (1-3) 
-oo 



Exact equivalence has been proved at least on a torus [2] . 
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Thus we can find Z = Y[n=o and the vacuum expectation value of ticf) may 
be calculated as 

11 N-l ^, 
( tr0'^) = -V (n|A^|n), (1.4) 

where A is the operator multiplying with A, and | n ) is the orthonormal basis defined 
by In) = exp (— I V{X)) Pn{X)/\^Ki. Henceforth the r.h.s. will be simply denoted 
by (A'^). 

Though the potential V{X) = YlkLo 9^ should be fixed at one of the critical 
points at the double-scaling limit, at finite N we are free to alter it from the critical 
point, and find that 

5 {Pul^n) - V Onm Pm/ VK., (1-5) 

where all upper half components of O are empty. It is convenient to work on the 
orthonormal basis, and define Qnm = {n \ \ \m). Qisa symmetric matrix, and 
8V{Q) =0 + *0. Thus O can be identified with [5V{Q)]-, where [ ]± denotes the 
projection such as \ [ ] diagonal + [ ] upper/lower triangular- Heuce the Variation of Q 
can be obtained as 

5q = [q\ mm^-wm-)]- a-e) 



Toda Hamiltonian flows defined by Adler [3] are simply 

[Q,[Q'=]+] = -[Q,[Q*]_], (1.7) 



dgk 

and from 5V{X) = ^(A — 5X) — V^(A), one also finds the discrete string equation 

[p,g] = i, (1.8) 

with P = li[V'{Q)]+ - [V'{Q)]-). As proved by Martinec [3], the Toda flows are 
commutative, and also preserve the string equation. 

In terms of the Ward identities under the Toda flows, Mironov and Morozov 
[5] have found the Virasoro constraints 

L„Z = 0, (n>-l) (1.9) 

where 

^1=1 dgi+n ^1=0 dgi dgn-i 

satisfies the Virasoro algebra [L„, Lm] = (n — m) Ln+m for n, m > —1. 



1.2 Double- Scaling Limit 



The GL(cxd) Toda hierarchy is transmuted into the KdV hierarchy at the 
double-scaling limit, and more generally, it has been shown by Douglas [6] that 
there appears the generalized KdV hierarchy in the one-cut hermitian multi-matrix 
models. The degree n of the orthogonal polynomials may be described in terms of a 
continuous parameter t at the double-scaling limit, such that n/ (3 — 1 — at* Then 
the pair Q and P, which are the matrices defined similarly for the first piece of the 
matrix chain, scale to the scalar Lax pair [7] as /3 — > oo and a — > simultaneously. 
In a critical multi-matrix model, for instance, one can find the Lax operator L = 
d1 -\- Uq-2 d1~^ -h . . . -h and (1.8) scales to the continuum string equation 

[(L^/'^)+,L] = l, (1.11) 

where ttq-2,- • • ,wo are functions of t called the scaling functions, and ( )-|- denotes 
the ordinary differential operator part of the pseudo differential operator L^/*?. Sim- 
ilarly to the Toda hierarchy, the generalized KdV hierarchy has infinite number of 
commutative Hamiltonian fiows 

^ = [(L^/«)+,L] = -[(LP/«)_,L], (1.12) 

and in the matrix model, these generalized KdV flows can be generated by the 
scaling operators ak{Oa) with p = kq + a + 1 such that 

^ = MOa))= dt{t\LP/'^\t), (1.13) 

where F = logZ, and the residue i?fc,a = (^| L^/'^ |t) is a differential polynomial 
of the scaling functions satisfying certain recursion relations [8-10]. Integrating 
the string equation (1.11), one also finds the Virasoro and the W„-constraints {n = 
3, . . . , g) on Z, and thereby Z can be identified with the r-function of the generalized 
KdV hierarchy [9-11]. Since the integrability appearing in the continuum is merely 
a consequence of the integrability in the discrete, it has been conjectured that the 
Virasoro and the W-constraints can be derived from the discrete Ward identities 
at any scaling limit. This has been proved at least for the one-cut hermitian one 
matrix models [11,12]. 

These constraints have been already found in the 2d topological gravity cou- 
pled with the ADE series of the topological minimal models [8,13], and very recently, 
the same Virasoro constraints has also arisen in the Kontsevich model [14,15]. 



Note that 1 — N/ (3 is the cosmological constant. 
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1.3 Overview 

The one-cut family of the hermitian matrix models has demonstrated their rich 
structures in the generalized KdV hierarchy and W-constraints. In this paper, we 
shall study the two-cut family of the hermitian one matrix models [16-18], and look 
for similar structures. The procedure developed in ref. [12], which is easier to handle 
and more straightforward than the one used in ref. [11], will be applied to the two- 
cut family; and consequently, the continuum Virasoro constraints will be derived 
from the discrete Ward identities (§3). The result agrees with refs. [17,18], and 
the meaning of the additional parameter a appearing in the continuum Virasoro 
constraints will be clarified. Eventually one can find the nonlinear Schrodinger 
(NLS) hierarchy for the 2x2 matrix Lax operators [19,20], and obtain the mKdV 
hierarchy as the reduction into the even potential models [17,18]. The invariance of 
the string equation under the NLS flows will be proved in §4. HoUowood et al. [18] 
have obtained what they call the Zakharov-Shabat (ZS) hierarchy by rotating the 
anti-hermitian model to the hermitian, and from the ZS hierarchy one can get both 
KdV and mKdV hierarchies by reduction. We will investigate this rotation further 
in §5.1 according to the mapping we found, and construct a rotation from the NLS 
hierarchy to the ZS hierarchy. Finally the quantization of a reported in ref. [18] 
will be interpreted as the analyticity of the model at the continuum limit e — > by 
means of the rescaling of e (§5.2). 

Generally one can show that the family of matrix models is governed by two 
kinds of distinct integrable structures in the discrete and the continuum, and the 
investigation on the relation between these two structures is not just of mathemat- 
ical interest, but is essential in order to understand the physical outcome from the 
integrability in the continuum. As we will examine later, in order to reach a gen- 
eral point (ti,t2, • • must perturb the matrix model by the associated flows, 
which may or may not be allowed within the original model. Especially the flow 
connecting two different critical points must deviate from the matrix model. This is 
because even though the critical points are individually realized by the critical ma- 
trix models, in between one cannot find the corresponding matrix models, in other 
words, if we impose the physical conditions on the solution of the string equation 
in order to extract the physics from the formal system like the generalized KdV, 
the NLS, or the ZS hierarchy, those solutions generally become unstable under the 
associated flows [21]. Unfortunately the mapping itself is not strong enough to solve 
this reduction problem, but it will provide a foundation to get further insight. 
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2. Two-cut Models 

2.1 Hierarchical Criticality 



We will determine the critical potentials of the hermitian one matrix models, 
and classify them into hierarchies. The string equation (1.8) governs the matrix 
model completely, but here we use an alternative, but equivalent pair of equations 



{n-l\V'{X)\n) = ^, 
(n| V'{X) \n) = 0, 

which we also call the string equation. In terms of the recursion relation 

AP,(A) = Pn+l{X)+SnPn{X) + RnPn-lW, 

one can show that 



X=^ ^xVh = z + S + R, 
Vh 



(2.1a) 
(2.16) 

(2.2) 
(2.3) 



where R \n) = i?„ |n), S \n) = Sn \n), h \n) = |n), z \n) = |n+l), and 
z^ \n) = \n — 1). Due to the hermiticity of A, (2.1a) may be rewritten as 



{n\l--XV'{X)\n) = l- 



n + 



(2.1c) 



which is more convenient for our present purpose. Now assume that Rn and Sn 
converge as n — > /3, and denote those limits by and 6, respectively (the positivity 
of Rn follows from Rn — hn/hn-i), then (2.1c) and (2.1b) may be evaluated by the 
contour integrals around z = 



r dz 

J Co ^ 



l-lx{z)V'{X{z)) 



0, 



/ ^V'{X{z)) = 0, 

JCo ^ 



(2.4) 



with X{z) ~ z + b + 0? jz. We can rewrite (2.4) further in terms of 2; = [A — 6 — 
v/(A- 6)^ -4a2]/2 such that 



\--XV\X) 



[ dX [(A - bf - 4a^] 
[ dX [(A - - 4aY^ V'{X) = 0, 



0, 



(2.5) 



where the contour Coo must be large enough to enclose all singularities on the 
complex plane. Noticing that 

/ dr]{r]^ -4a^)-^ r]"" = 0, (2.6) 
for m odd integral or negative even integral, it is easy to obtain a general solution 



1 

2- — ^2.7) 

F{X) = c(A) (A - A+)^ (A - A_)i 



l__Ay'(A) = [AF(A)]+, 



Here \± — b ± 2a, and [ ]_|_ denotes the principal part (including the A^ term) of 
the Laurent expansion about A = oo. c(A) is an arbitrary entire function satisfying 
ResA=oo F{\) = —1 in order that the A° term of AF(A) equals 1. Note that the first 
line in (2.7) may be expressed as V'{X) — —2 [F(A)]^. F{X) controls the eigenvalue 
distribution of A through (see e.g. [22]) 

p(A) = - lim ImF(A-ie). (2.8) 

TT 6^0+ 

Hence c(A) must be negative along the support of p(A), which turns out to be the 
interval (A_, A_|_) separated by the even order zeros of c(A) (no odd order zeros are 
allowed along (A_, X+)). One might think c(A) could have cuts along (A_, A-|-), but 
this cannot happen because of (2.5) and the negativity of c(A). We will see, for the 
critical potential V^(A) given by (2.7), that (2.1a) and (2.1b) scale to the continuum 
string equation at the double-scaling limit, and this scaling behavior is completely 
governed by the eigenvalue distribution at the edge, i.e. the zeros of F{X). 

Suppose c(A) has a zero only at A = A-|-, then c(A) may be expanded as 

where must satisfy ^fc>i (20)^^^^ ''(fc+2)!" '^k — —1- This is the family of one-cut 
models, in which the m*'* critical potentials can be obtained by choosing < 
and Cfc = for A; < m, where the negativity of follows from the negativity of 
c(A). More generally, c(A) may have zeros at both ends such as 

c(A) = E,^,fe_>o - ^)'^ (^ - ^-)'"' (2.10) 

which is the doubling studied in ref. [16]; in particular, k+ — k- gives the family 
of even one-cut models. All of these one-cut families yield the KdV hierarchy at 
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the double-scaling limit, whereas the niKdV hierarchy and more generally the NLS 
hierarchy may be obtained from the two-cut family as follows. 

Suppose c(A) has an even order zero at A = 6 

^(^) = E,>l^2fc(A-6)2^ (2.11) 

then C2k niust satisfy X^fe>i (2a^)'^"'"^ ^1^+1)1' = —1- These are the two-cut mod- 
els studied in ref. [18], in which the 2m*'* critical potentials* are given by C2m < 
and C2fc = for A; < m; in particular we caU C2k = - (2a^)iTi\2m-i)!! ^k,m the 2m*'* 
critical point. More general critical model can be obtained by taking the product 
of these two kinds of critical models in the same manner as the (/c_|_, /c_) model has 
been constructed from the doubling of the A;+*'* and A;_*'* critical one-cut models in 
ref. [16]. 



2.2 String Equation 

We will fix our notation for the two-cut models by following Crnkovic and 
Moore [17], and give a heuristic argument to derive the continuum string equations. 
The 2m*'* double-scaling limit may be defined by e — > in 

i?„ = l + (-l)-26/(t) + ..., 

Sn = b+{-ir2eg{t) + ..., ^ ■ ^ 

with X — n/ (3 — 1 — e^'^t, and k"^ (3 e^'^'^^ = 1 (a has been normalized to 1). For 
the orthogonal two component vector defined by 



(2.13) 



P+(A) = exp (-^ V(X)) i-iy P2i{X)/Vh^i, 
P27+i(A) = exp (-| V{X)) {-ly P2i+i{X)/Vh^i, 
A may be represented as a 2 x 2 matrix operator 

At the scaling limit, -P^(A) scales to K,^/e{t\X)^ with | A ) being a two-component 
vector, while both A and A scale as 

(t| A |A) = (t| 6E + 2eL |A). (2.15) 



No odd critical potentials are allowed within the hermitian models. 
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Hence we will not distinguish A from A. E is the identity 2x2 matrix, and L is the 
Lax operator L = —ia2 dt — (Jif + (T^g with ct^'s being the Pauli matrices. 

For the critical potential given by (2.11), the leading part of (2.1a) is trivially 
satisfied due to (2.4), and the rest of (2.1a) and (2.1b) can be combined into the 
string equation 

{t\V'{X)\t) = -ie^"'t{t\t), (2.16) 

where the factor —i is necessary to compensate for the extra i factor contained in 
{t\t) (see §2.3). Based on the same argument as ref. [12], all {t \ X"' \t) with 
n negative vanish. Basically this is because of the analytic continuation from n 
positive to negative. Thus we can remove the restriction [ ]+ from V'{X), and find 
that 

t{t\ \t)+A'^+^C2m{t\ L^"^ 10 = 0- (2-17) 

This holds for any 2m*'* critical potential, since the higher order terms disappear 
as e — > 0. 



2.3 Diagonal Kernel of the Resolvent 



The explicit form of the residues for the general matrix Lax operator L can 
be obtained from the resolvent L — XE, which has been studied in ref. [20] and 
their theorem 6 provides us with the following results. The Laurent expansion of 
the diagonal kernel of the resolvent about A = oo starts from the zero-th order, and 
for k > —1 



{t\L''\t) = - {Fkas + Gkai + HkE), (2.18) 



where 



-^fc+i - gHk- - Gfc, 

Gk+i = -fHk + ^F;^, (2.19) 

with ' denoting D = df. First two levels can be determined as F_i = G-i = 0, 
H-i — 1, Fq — g, Go — — /, Hq = by assuming that / and g are t independent 
and directly evaluating the diagonal kernel as 



^ ' L - A S ' ' e-^o+ J_^ 2t: \ -f + ip -g-X + ie 



(2.20) 



8 



Higher levels may be obtained from (2.19), where the first few are 

F^ = -\9" + \ 9if + 9'), G, = \ f" - \ fif +g'), H, = \ (gf - fg'), 
Fs = -l r + jf'if + 9% Gs = -l g'" + jg'if + g% 



2.4 Reduction to the mKdV Hierarchy 

For an even potential model, one finds 5'^ = and therefore g = 0, for 
which (2.19) yields F2k+i = | S'^^ifi G2k = -Sk[fi H2k+i = Rk+i[u] - \ S'^[f] = 
R}~j^i[u\ + \S'iJ[f]^ and otherwise zero, where u and u are defined by the Miura 
mapping u = p + f and u = p — and Sk[f] and Rk[u] are the k^^ Gel'fand- 
Dikii differential polynomial in the mKdV and the KdV hierarchies, respectively. 
The definitions are 

Sk[f] = [(/ - ^ D) U+\d) Df f^[-\D' + fD-'f D]' /, 

Rk[u] = [D-'{f+lD)D{f-lD)]'^l = [-\D' + lu+^-D-'uD]H. 

(2.22) 

The 2m*'* string equations in the NLS hierarchy are therefore 

tg + A'^+^C2mF2m = 0, (2.23a) 
tf-4^+'c2mG2m = 0, (2.236) 

with C2m < 0. (2.17) also provides the third string equation H2m = 0, but since 
^2m ~ follows from (2.23), and H2m is odd under either / — > — / oi g ^ —g 
while any of (±/, i^f) satisfies (2.23), the third string equation may be obtained by 
integrating (2.23).* 

We will conclude this section with a few comments. Since the scaling behavior 
does not depend on the parameter 6, wc omit h hereafter, while keeping g general. In 
the sphere limit k = 0, or as t — oo, all t derivative disappear, so that (2.20) provides 
exact relations, such as = g {p+g^r, G2m = -^^^ / {^ + 9^^- 

The string equation at the 2m*'* critical point is thus given by 

t = 2{m + l){f+g^r. (2.24) 

* The integration constant does not vanish if the odd modes are mixed in. 
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For an even potential, the system is reduced by = 0, and only (2.23b) remains 
nontrivial, which may be written down as 



at the 2m*^ critical point. The same string equation has been found also in the one- 
cut family of the unitary matrix models in ref. [23]; furthermore it is proved that 
(2.25) has a unique, real, pole-free, solution consistent with the above asymptotic 
behavior (2.24) for t ^ oo [24]. 



3. Discrete and Continuum 

3.1 Mapping from Toda to NLS 

The partition function (1.1) is a function of coupling constants through the 
potential V{z) = J2k^=o 9k i and satisfies the differential equation 



J{z)Z 



\v'{z)-{W{z)) 



Z. 



(3.1) 



Here {W{z)) = {(3 ^ ir[{z — (j)) ^]) is the generating function, and J{z) = 



,-k-l 



Jk is defined by Jk = P d/dgk for A; > 0, and J-k = ^ gk for k > 1 [12], so 



ez 



that Jfc satisfies 



[Jk, Ji,] = -^P '^k6k+i,o. 



(3.2) 



Re-expanding J{z) as J{z) = J^fcez ^ ^ {z'^ — ^) ^ Jk, one can find that 

dz 



J2k — 



J2k+i= I ^{z'-4)-iz''J°''%z), 

J Coo 



dz 



(3.3) 



where J{z) is divided into odd and even functions of z, denoted as J'^'"^'^{z) and 
J°'^'^{z), respectively. Note that the superscripts of J are named according to the 
even/odd nature of the potential, not its derivative. J°'^'^{z) disappears from 
because of the cancellation along Cqo, and so does J^'"^'^{z) from J2k- From (3.2) 
and (3.3), the commutation relations for Jk can be calculated as 



Jk , Jl 



= P- 



-k5k+i,o -2{k-l) 5k+i,2 



(3.4) 
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According to (2.11), V' may be expand such that 



(3.5) 



where the 2m*'* critical point is defined by Ck = ~ 2"'+U2m-iy.\ ^k,2m- From (3.4), 
the non-positive modes J/.; k = 0, —1, —2, . . . commute with each other and contain 
no differential operators except for d/dgo in Jq, which may be replaced by —x in 
terms of I3~'^ d {log Z)/d go = —x. Therefore these may be considered as a new 
coordinate system for the space of the coupling constants instead of c^'s or g^s. 
These are related to Ck such that 



def 



J_fc 2{k + l)pk = 4:Ck- Cfe-2, 
where c_2 = x — 1, since Res^=< 



{k > 0) 



(3.6) 



Efc>-i Cfcz'=(z^ -4)^] = -1. On the other 
hand, the positive modes Jk, k = 1, 2, . . ., are the pure differential operators w.r.t. 
Pk, namely 



Jk = P 
1 



-2 



k 



d 



d 



8 dpi 



4{k+l)dpk dpk-2 
d 



{k > 2) 



(3.7) 



This is not the case if we write down J^'s as differential operators w.r.t. ^^'s, i.e. 



{k>l) 



where 



JkZ = - 



-1 / \k-l 



(3.8) 



(3.9a) 



and 



_ 1 /■ dz 



4)2 ^'=-iy'(^), 



(3.96) 



linear combination of gi-, 



Since J J, is a t independent number, we can absorb into the integration constant 
of the trace ( ) for each J^, and obtain 



JkZ ^ -(3-^ {X^-^ (X^ - 4)^) Z. 



(3.10) 



Defining the scaling parameters by 
4c_i = 2e^"*+^ 



9, 



4co -c_2 = e^^to, 



(3.11) 



4 Cfc = e 



2m— k n—k 



2-'=(/c + l)tfc, 
11 



(fc>l) 



the 2m*^ critical point of the two-cut model is given hj q = 0, to = t and tk = 
- ''"^(Tj^iti^' ^k,2m- As e ^ 0, Cfc_2; > 1 in (3.6) and d/dpk; k >2m (3.7) 
become negligible due to the relative e factor, and hence scales as 

Jfc = -K4e2m+/c2fc-i^^ (;.>2) 

Otk-2 

J_fe = e2"^-'=2-'=(A;+l)tfe, (k > 0) 

where from (3.9) one can obtain 



2^v- 2'+i(2Z-l)!! 



K 



2 



d _i d 



(3.13) 



Using (3.12) in (3.10), one can find the commutative fiows in the NLS hierarchy 
such that 

^2 1^ ^ ^-k-i ^-k-i0k+i (p - 4)^), 
atk 

/oo 
dttT{t\L^+'^ \t), (3.14) 

/oo 
dtHk+ii 

for /c > — 1, where F = log Z and the trace on the 2x2 matrix is necessary in order 
to sum up the contributions from even and odd polynomials. Since dF/dt^i = 
follows from Hq = 0, the t_i derivative vanishes in Ji. Note that a appears in 
these relations because it cannot be absorbed into the integration constant, a will 
be identified with the integration constant of the third string equation, and proved 
to be independent of t^s in §4.2. 
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3.2 Virasoro Constraints 



Multiplying (3.1) once again with J{z), one can find the discrete Ward identity 



\12] 



■.J{zf: Z^-A{z)Z. 



(3.15) 



where : : denotes the normal ordering, and A (2;) is an entire function defined by 



A(.)=(v'wf_4((i:(fi5p)|). 



(3.16) 



Hence Z satisfies 



dz 
27H 



: J{zf : Z = 0, (n > -1) 



(3.17) 



for which one can find the discrete Virasoro constraints (1.10) by calculating the 
contour integral. Alternatively (3.17) may be rewritten as 



dz 



I ^ : Jizf : Z = 0, 
f ^ ^ : j(^)2 : z = 0, 

I ^z-+'{z'-A):J{zr:Z = 0, (n > -1) 
and calculating the contour integral we can replace (1.10) with 



E 



fcez 



(E. 



-21- 



Z = 0, 



ez 



. J — fc J'} 



k 'Jn+k-2 



:) Z = 0. 



From (3.3) and (2.6), we can get 



Coa —Cq 

dz 
27H 



— ^(^2-4)-^ J"(^) Z, 
2ni 



{z'^-A)-^ 



1 



zV'{z) - 1 



(3.18a) 
(3.186) 
(3.18c) 



(3.19a) 
(3.196) 
(3.19c) 



(3.20a) 
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/ (3.206) 



where the r.h.s. of both (3.20a) and (3.20b) vanish automaticaUy, since the critical 
potential (3.5) always satisfies (2.5). Hence (3.19a) and (3.19b) are satisfied trivially 
before taking the scaling limit. In terms of the scaling of in the vicinity of the 
2m*'* critical point, (3.19c) scales to the continuum Virasoro constraints 

L„Z = 0, (n>-l) (3.21) 



d 

— 2 a t, 



where 




L-i 


^ \00 


Lo 








Ln 





dtk-i 

9 -4 2 

dtk 
d d_ 

dtk+i dV 

+ a 



(3.22) 



fe+n 4 -^fc=l dtn-k-1 dtn-1 

These constraints have been found already in refs. [17,18] by integrating the string 
equations 

J]^^^(fc + l)t,Ffc = 0, 
T.k>o if' + ^)i>^Gk = 0. 
Finally, defining (afc) = J* dt Hk+i, we arrive at 



(3.23) 



, {k + l)tk{(Jk-i) - K-^at = 0, 

k=l 

2 

■r — ^00 „ Ct 

J2k=o (^ + 1)*'^ — = 0, 

(A; + 1) tfc (afc+i) + a (ao) = 0, (3.24) 

k=0 

/ l)tk{crk+n) + o:{<7n-i) 

/ V, . (("^fc-l CTn-fc-l) + (o-fc-l) (c^n-fc-l)) = 0. 



14 



4. NLS Hierarchy 

4.1 Flows and String Equations 

Gcl'fand and Dikii [20] have obtained the Lax pair for the general matrix 
differential operator L in terms of the diagonal kernel R{z) = {t \ {L — z E)~^ \ t). 
For our first order 2x2 matrix operator L, their theorem 7 can be read as 

y Pk z-^-^ = R{z) a2{L-zE)-^ = {L-z E)'^ a-, R{z), (4.1) 

where the last equality follows from their theorem 6. First few of are 



Pi = \a2{D'-\ {D, fas + gai} - ^ {fas + gaif], (4.2) 
P^ = l[D'-l {^'' M + gcTi} - I {D, {fas + ga,f}], 



and the k*'^ equation in the NLS hierarchy may be defined as 

2[L,Pk+i]=2Gk+ias-2Fk+iai. (4.3) 



dL_ 
dtk 



Here the k*'' NLS fiows 



dl 
dtk 
dg_ 
dtk 



2 -Pfc+i, 

(4.4) 

2 Gfc+i) 



are all commutative for k> —1. 



The original definition in ref. [20] has no flows for k = —1 due to the extra 
{k + 1) factor appearing in the r.h.s. of (4.4), which we have been normalized to 
the definition of tk- Both are consistent with the previous definition based on the 
mapping from the Toda flows. In fact, by taking the t derivative twice on (3.14), 
one can flnd that 



^^^^Kp^ = 2 i^Ui = 4 (/ Fk+, + g Gk+i), (4.5) 
dtk 

for k > —1, which agrees with (4.4), while the (— l)*'^ flow vanishes on + g^ 
anyway. Actually these two deflnitions are equivalent due to the factorization of 
the t-i dependence. This can be done as follows. First introduce the complex 
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notation ifj = f + i g and Uk = Fk + i Gk {k > —1); thus the t-i dependence of '4> 
can be factorized as '0 = g-Sjt-i |-,y solving the (—1)*'^ equation 

= 2^, 



(4.6) 

Rewriting (2.19) with the complex notation* 

Uk+i = -ii^Hk+'-Ui, ^^^^ 

and using the induction w.r.t. k starting with U-i = —iij), one can find that all 
Uk have the same t_i dependence as i/;, namely Uu = e~^**-i iUk)o, and t-i in- 
dependent '00 and (C/fc)o again satisfy the same recursion relations as (4.7), where 
H^s are t-\ independent by virtue of (4.5). Therefore we can conclude that the 
t-i dependence is redundant, and has a global phase ambiguity. We will leave 
the (-1)*'* flow non-vanishing for convenience. 

It is noteworthy that the string equations (3.23) can be rewritten as 

[M,L]=E, (4.8) 
by defining M = X]fc>o + -'-) ^k- This is the scaling from (1.8) instead of (2.1). 

4.2 Invariance of the String Equations 

We will prove the invariance of the string equations under the NLS fiow. 
Defining 

the string equations (3.23) may be expressed as 

K_i/ = K_i^ = 0, (4.10) 

and integrating K_i (/^ + g"^) = once w.r.t. t, one can obtain the third string 
equation 



* The first equation in the NLS hierarchy is — —t\)"j1 + the nonlinear 

Schrodinger equation. 
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where the integration constant c may not vanish, and by taking the t derivative 
once on L_i Z = 0, one can identify it as c = a. 

The string equations (4.10) and (4.11) are invariant under the NLS flows if 

(Z + l)Fi + K_iFi+i =0, (4.12a) 
{l + l)Gi + K_^Gi+i=0, (4.126) 
{I + 1) Hi + Hi+i = 0, (4.12c) 

hold for I > —1. Here we have used the tk independence of a to get (4.12c) 
from (4.12a) and (4.12b). da/dt-i — da/dto — foUows directly from the string 
equations, and da/dtk = {k : even) follows from the definition of a, while 
da/dtk — {k : odd) is a result of the fine-tuning of the odd couplings to make 
a finite. Now we can prove (4.12) by the induction w.r.t. / starting from / = — 1, 
which is nothing other than the string equations (4.10). Suppose that (4.12) holds 
for Z < A; — 1, then the recursion relations yield 

i^-i =9K_,H,-l (K-^ G,y + ^ ^ G„ 
= -(/c + l)Ffc. 

(4.12b) may be proved similarly. 

In the presence of (4.11), we can derive alternative string equations by taking 
the t derivative once on (4.10); namely 

Kof = -f, Kog=-g, (4.14) 

where 

Taking the t derivative n times on (4.14), one can find that 

= _(n + 1)/H, 

(4.16) 

i^o^("^ = -(n+l)^^'^>. 

It is convenient to define "dimensions" by d{f) — d{g) = — 1 and d{th) — k+1. Since 
Kq is a first order differential operator, i.e. Kq {AB) = {KqA)B + A{KqB), is 
therefore an operator counting the dimensions of the differential polynomial A[f^ g] 
if / and g are the solution of (4.14); namely 

KQA[f,g] = d{A)A[f,g], (4.17) 

where A[f.,g] may have explicit tk dependences except for t-i. It is then easy to 
prove the invariance of (4.14) under the NLS flows directly. 



17 



By restricting the potential to be an even function of A (6 7^ 0), we can find 
a constraint g = 0, and consequently the NLS hierarchy is reduced to the niKdV 
hierarchy with t2k-i = for /c > 0, where the odd flows are all frozen because they 
vary the constraint ^ = 0. Only the second string equation K_i g = remains 
nontrivial in this reduction, and one can obtain the mKdV flows and the string 
equation 



dl_ 
dt2k 



-s'M 



{k > 0) 



i2k+l)t2kSk[f]=0, 



(4.18) 
(4.19) 



where the string susceptibility is given by df F 



K 



-2 f2 



f . It is more appropriate 



to describe (4.19) as the once integrated form of Kq f = — /, which is manifestly 
invariant under the mKdV flows. 

5. ZS Hierarchy 

5.1 Rotation from NLS to ZS 

One can deal with the anti-hermitian matrix model in the same manner as 
we did with the hermitian matrix model in §2. HoUowood et al. [18] have already 
studied the anti-hermitian matrix model with real couplings, and found what they 
call the Zakharov-Shabat (ZS) hierarchy. Here we will re-investigate this model in 
terms of the mapping we have found. By rotating the anti-hermitian matrix (j) to 
the hermitian matrix (f) = i(f), one can flnd that the hermitian matrix model has 
a complex potential, namely gk = {—i)^gk) where ^^'s are the couplings of the 
anti-hermitian model, and all real valued. Since the eigenvalues of (j) distributes 
along the imaginary axis, the limit values of Rn and 5"^ also rotate as a = — z a and 
h = —ib, where a and b are the limits deflned previously in the hermitian model, 
and consequently the scaling functions are given by / = — / and g = —ig with / 
and g being real valued. Defining the coupling constant of the hermitian matrix 
model by Ck = (— i)^"'"'^ (cfe e R), and the spectral parameter by z = i z (2; e R), 
one can find that 



k>-l 



Ck 



{z - bf+^ {z - X+)i {z-X-) 



2 tz 



(5.1) 



in which the k^^ critical potential is either real or pure imaginary depending on k 
even or odd as we expect. C2k must satisfy the constraint J2k>i 0,'^)'^'^^ '^'{k+i)\' ^2fe 
= —1, but unlike the previous model, C2k-i does not have to be zero, because the 
pure imaginary potential does not contribute to the eigenvalue density. Now the 
Lax operator 

L = —i L = —i a2 D + ai f — i g , (5-2) 
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turns out to be complex, thus we need an SU (2) rotation in order to make it pure 
imaginary, i.e. 

L ^ a^D + ia2 f + aig. (5.3) 

Under the rotation from to (— i)'^"'"^ c^, the NLS hierarchy is rotated to the 
ZS hierarchy by tk tfc, where we have changed tk to (—1)'^"'"^ tk for convenience. 

Note that to = and the 2m*^ critical point is given by tk = ^~^(2m+i)!r^''' ^k,2m- 
Defining new differential polynomials in such a way that the fiow equations appear 
to be the same, namely 

Fk ^ {-if^^ Fk, Gk ^ i-i)'' Gk, Hk ^ (-1)''+' Hk, (5.4) 
with F-i — G-i — 0, H-i — 1, Fq — g, Gq — f, Hq — 0; we can find the residue 

(t|L'=|t) = ^ (Ffeai+iGfeda + ^fe^;), (5.5) 
with the recursion relations 

Fk+i = gHk + G'ki 

Gk+i^fHk + \Fi (5.6) 
^-H'k=gGk-fFk. 

Now all negative signs in the Virasoro constraint (3.22) change to the positive signs 
because of the factor, but (3.24) does not change because 

«2-^ = 2 / dtHk+i. (5.7) 
Otk Jt 

In terms oi if) = f — g., if} — f -\- g., Uk = Fk — Gk, and IJk — Fk + Gk, the ZS fiows 
and the string equations may be rewritten as 

^ = 2Uk+i, ^^2Uk+i, (5.8) 
otk otk 



E,>o + 1) = 0' Efc>o + 1) 

V {k+l)ikHk = c. 



(5.9) 



For the even potential, again g = reduces the ZS hierarchy to the mKdV 
hierarchy, for which the differential polynomials are F2k+i — (— l)'^ h ^k[f]-i ^2k = 
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(-1)'= Sk[f], H2k+i = (-1)^ (I S'j^if] - Rk+i[P + f']), and zero otherwise, hence we 
can find the same mKdV flows as (4.18) and the same string equation as (4.19) by 
assigning t2k = (~1)'^ ^2fc7 and also the string susceptibihty is again d^F = p . 
Note that the 2m^^ critical point, t}~ = ~ ^ {2ra+r)y^}^ ^fc,2m; is precisely same as 
before. 

HoUowood et al. [18] have found an alternative way to freeze odd flows, namely 
t/J = e^*-i and u = '^'0, for which one can find that 

U2k = (-1)' (^ - u) Rk[u], U2k = H2k-i = {-!)'' Rk[u], 

U2k+i =2H2k+2 = {-l)''R'k+iH ^^^^^ = °' 

with Rk [u] being the k^^ Gel'fand-Dikii differential polynomials of the KdV hierar- 
chy, therefore this is the reduction into the KdV hierarchy 

=2R'^^M. (5.11) 



dt2k 

V (2/c+l)t2fei?feM = 0, (5.12) 

'fc>0 

where the string susceptibility is given hy F = u. This result agrees with 
the one-cut family of the hermitian models (renormalization of t2k is necessary to 
adjust the critical points). The KdV string equation is the once integrated form of 
the third string equation with c = |, and is therefore manifestly invariant under 
the KdV flows. 



5.2 Quantization of c 

The NLS hierarchy may be obtained from the two component KP hierarchy 
by the so-called reduction procedure, or alternatively, from the hierarchy of soliton 
equations associated with the A^^^ Kac-Moody algebra, for which two different ver- 
tex operator realizations (the principal and the homogeneous pictures) respectively 
lead to the KdV and the NLS hierarchies [25]. In the homogeneous picture, the 
vertex operator has one free parameter 7 due to the zero mode of the scalar field, 
where 7 could be any element of the finite root lattice, and consequently the tau- 
function depends also on 7. For A^^\ one finds 7 = mcti (m G Z), and the scaling 
functions '4> and t/S may be assigned to be 

^ = «2+4cIm+l^ ^ = ^2-4cZm-l^ V = log T^, (5.13) 

for which one can get the NLS and the ZS hierarchies by taking particular real 
slices, i.e. ip = ip* (resp. '0, -0 e R). 
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By following this tau-function formalism, HoUowood et al. [18] have found 
that the Virasoro constraints L^(c) r^, = (n > — 1) produce 

L„(ct1)t^±i = 0, (n>-l) (5.14) 

which strongly suggests the quantization of c. 

Actually c must be quantized as c G Z/2 due to the following mechanism. If 
we rescale e to ae, while keeping the model and k fixed, tk changes to a^"*"^ tk- This 
can be seen most easily in the first line of (3.14), where F and ( ) do not change by 
assumption. Since (4.14) holds as well in the ZS hierarchy, and the t_i dependence 
may be factorized as ip = e~^*-^ Vo and ip = e^*-^ V^o, (and ip) has the anomalous 
dimension —2c (resp. 2c), i.e. 

E.>o(^+l)'~^f^ = -(^ + 2c)V^, 

fj (5.15) 

Therefore if 2c is not integral, the continuum limit e — > becomes non-analytic 
with respect to e, and hence c must be half-integral. 

(5.15) also explains the k factors appearing in (5.13), because the rescaling 
of to may be realized by the rescaling of k^. If the way Rn and Sn approach their 
limit values is scale invariant, c must be zero, which is the case for the reduction to 
the mKdV hierarchy. On the other hand, the reduction to the KdV hierarchy has 
been achieved by the constraint ^p — e^*-% thus the scaling dimension of ip must be 
zero, and hence (5.15) gives c = 1/2 which agrees with the previous calculation. 

We have derived (4.14) from the string equations (4.10) and (4.11), therefore 
any solution of the string equations (5.9) satisfies (5.15) for arbitrary value of c. 
However, there is no guarantee that (5.15) with arbitrary c is compatible with the 
constraint we require such as ffc = {k : odd). Above examples show that the 
compatibility actually breaks down. But of course, not all constraints quantize c, 
for instance, in the topological phase emerging at ti and = for k >2, the 
string equations have exact solution such as 

^ = -g-tV^"'exp(-2t_i + t-), V = 5^tr^^'exp(2t_i - ^), (5.16) 

with r + s = 1, which satisfy (5.15) and the third string equation for arbitrary c. 
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6. Discussions and Conclusion 



6.1 Hermitian Models and the NLS Hierarchy 

The 2m*'^ mapping we have obtained in the hermitian model is 

4co-c_2 = e2"^to, {to = t) ^^-^^ 

4cfc = e2^-'=2-'^(/c + l)tfc, (A;>1) 

which maps the 2m*'* critical potentials, i.e. C2m < and Cfe = {k < 2m), to the 
subspace t2m < and tk = {k > 2m). At e = 0, only the 2m*'* critical point 
is admissible, because c/- = {k < 2m) and = {k > 2m) respectively imply 

tfe = {k < 2m) and Cjt = (/c > 2m) at e = 0, thus t2m = (2m+i)!! fro™ ^^e 

constraint. Owing to the positivity of the density of the eigenvalue distribution, 
must vanish for k odd, so that no odd critical points are allowed within the hermitian 
model; nevertheless we may turn on a, ti, . . . , t2m-i: since c_i = . . . = C2m-i = 
hold as e — > 0. To this extent, we can perturb the system away from the 2m*'* critical 
point and generalize the system to allow C2m > as well as the odd critical points. 
(6.1) should be understood in this sense, and also we extend it to the (2m — 1)*'* 
mapping, where any point of the subspace t2m-i 7^ and = (/c 7^ 0, 2m — 1) 
may be chosen as the (2m — 1)*'* critical point. 

Among these "time" parameters, the (—1)*'* parameters a and t-i are very 
special, a may be identified with the integration constant of the third string equa- 
tion, and in order to get finite a we must fine-tune the odd couplings such that 

5^fe>i ^ <^2fe-i = c-i, or equivalently, 1a + Y.k>i i^-^ik+iy. *2fe-i = 

however we do not have any flows associated with a, instead the ( — 1)*'* flow changes 
the parameter and its numerical value is independent of any of the coupling 
constants. The reason of the quantization of c = a is therefore not because of 
the momentum zero-mode quantization such as 2c = —K^d/dt-i., but because of 
the analyticity of the continuum limit discussed in §5.2. This is somehow analogous 
to the quantization of the central extension c in the Kac-Moody algebra, where c 
commutes with its dual d (the zero mode of the Virasoro algebra). 

The image of the generalized m}^ critical potentials consists of two discon- 
nected components depending on tm negative or positive, which we call the negative 
(positive) side of the m*'* leaf. The NLS flows from the flrst to the m*'* order deflne 
smooth coordinates on each side of the m*'* leaf, where the even critical point is 
located on the negative side of the even leaf, so that the odd leaves and the positive 
side of even leaves are unreachable by the matrix model. The 2m*'* flow cannot 
pass through t2m — due to the constraint we have, and because of that, if we im- 
pose the "physical" conditions to the solution of the string equations, i.e. real and 
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pole- free along the positive real axis and proper asymptotic behavior for t ^ oo, 
that solution has an instability at t2m — 0. 



6.2 Anti-Hermitian Models and the ZS Hierarchy 



For the anti-hermitian model, (6.1) may be modified as 
2 c_i - 2^^^^ —--J C21-1 - e 2 a, 



4co-c_2 = e-to, {io = t) ^^'^^ 



45k = e^-^ 2-'^ (fc + 1) 4. {k > 1) 

Rotating NLS to ZS, we can find that the above arguments hold perfectly good, 
except for the fact that the anti-hermitian model can realize the odd critical points 
as well, so the odd leaves are reachable by the anti-hermitian matrix models. Con- 
sequently, at the (2m — 1)*'* critical point, Cfe < for the first non- vanishing even 
coupling constant and k must be equal to or higher than 2m. Furthermore it is not 
possible to have only one odd coupling constant C2m-i non- vanishing in order to 
get finite a. 



6.3 Conclusion 



We will conclude with a few comments. So far we have studied the mappings 
individually, but since these leaves cover the entire t-space, and also the structure of 
the NLS (ZS) fiows is universal, we may sew them and get a universal mapping. This 
universal system is however highly mathematical unless one imposes the physical 
conditions, and among the points of the t-space, only the 2m*'* critical point has a 
definite scaling dimension 7str = — 1/m. In fact, by solving the 2m*'* string equation 
asymptotically, one can find that 

^(*) = E,>o(^'^'"'"')'"'^'*' (6.3) 

which agrees with the scaling behavior of the continuum gravity coupled with con- 
formal matter. 

The quantization of c is necessary to make the continuum limit analytic, and 
it can be realized by fixing the odd flows by ^ = 0, or = e^*-^ with c = 
(resp. c = I). Combining these results with (5.14), it seems that there are no other 
inequivalent fixings to realize the half-integer c, but this is an open problem, as is 
the interpretation of the quantization of c topologically or field theoretically in the 
continuum quantum gravity. 

The first critical point in the ZS hierarchy has been studied in ref. [17] as the 
topological phase, but so far no topological models have been found to be equivalent 
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to this phase. Nonetheless this does not discourage us to search for new topological 
gravity, which may usher in a new geometry such as the topological phase of the one- 
cut hermitian model did in the intersection theory on the moduli space of Riemann 
surfaces [14,15]. 



Acknowledgements 

I would like to thank Prof. L. N. Chang for careful reading of this manuscript and 
stimulating discussions. 



References 

1. D.J. Gross and A.A. Migdal, Phys. Rev. Lett. 64 (1990) 127; Nucl. Phys. 
B340 (1990) 333; E. Brezin and V.A. Kazakov, Phys. Lett. B236 (1990) 144; 
M.R. Douglas and S.H. Shenker, Nucl. Phys. B335 (1990) 635. 

2. M. Bershadsky and LR. Klebanov, Phys. Rev. Lett. 65 (1990) 3088. 

3. M. Adler, Inv. Math. 50 (1979) 219; E. Martinec, Commun. Math. Phys. 138 
(1991) 437. 

4. H. Neuberger, "Matrix Models and KdV," preprint RU-90-58 (1990). 

5. A. Mironov and A. Morozov, Phys. Lett. B252 (1990) 47; H. Itoyama and Y. 
Matsuo, Phys. Lett. B255 (1991) 202. 

6. M. Douglas, Phys. Lett. B238 (1990) 176. 

7. V.G. Dringel'd and V.V. Sokolov, J. Sov. Math. 30 (1985) 1975. 

8. E. Witten, Nucl. Phys. B340 (1990) 281; R. Dijkgraaf and E. Witten, Nucl. 
Phys. B342 (1990) 486; R. Dijkgraaf, H. Verlinde and E. Verlinde, Nucl. 
Phys. B348 (1991) 457; "Notes on Topological String Theory and 2D Quan- 
tum Gravity," preprint PUPT-1217 (1990), to appear in Proc. of the Cargese 
Workshop. 

9. R. Dijkgraaf, H. Verhnde and E. VerUnde, Nucl .Phys. B348 (1991) 435. 

10. J. Goeree, Nucl. Phys. B358 (1991) 737. 

11. M. Fukuma, H. Kawai and R. Nakayama, Int. J. Mod. Phys. A6 (1991) 1385; 
"Infinite Dimensional Grassmanian Structure of Two-Dimensional Quantuam 
Gravity," preprint UT-572-TOKYO (1990). 



24 



12. W. Ogura, Mod. Phys. Lett. A6 (1991) 811-818. 

13. K. Li, Nucl. Phys. B354 (1991) 711, ibid. 725. 

14. M.L. Kontsevich, Funk. Anal. Appl. 25 (1991) 123. 

15. E. Witten, "On the Kontsevich Model and Other Models of Two Dimensional 
Gravity," preprint IASSNS-HEP-91/24 (1991); S. Kharchev, A. Marshakov, A. 
Mirnov, A. Morozov and A. Zabrodin, "Unification of All String Models with 
c < 1," preprint ITEP-M-8/91 (1991); D. Gross and M.J. Newman, "Unitary 
and Hermitian Matrices in an External Field IL The Kontsevich Model and 
Continuum Virasoro Constraints," preprint PUPT-1282 (1991). 

16. P. Petropoulos, Phys. Lett. B261 (1991) 402; S. Dalley, "Critical Conditions 
for Matrix Models of String Theory," preprint SHEP-90/91-6 (1990). 

17. C. Crnkovic and G. Moore, Phys. Lett. B257 (1991) 322; C. Crnkovic, M. 
Douglas and G. Moore, "Loop equations and the Topological Phase of Multi- 
cut Matrix Models," preprint YCTP-P25-91 (1991). 

18. C. Nappi, Mod. Phys. Lett. A5 (1990) 2773; T. HoUowood, L. Miramontes, A. 
Pasquinucci and C. Nappi, "Hermitian vs. Anti-Hermitian 1-Matrix Models 
and Their Hierarchies," preprint IASSNS-HEP-91/59. 

19. V.E. Zakharov and A.B. Shabat, Func. Anal. Appl. 8 (1974) 226; 
LM. Krichever, Func. Anal. Appl. 11 (1977) 12. 

20. LM. Gel'fand and L.A. Dikii, Func. Anal. Appl. 11 (1977) 93. 

21. For instance, it has been shown that the physical solution is unstable under the 
flow from the second critical point to the pure gravity (the first critical point) 
in the one-cut family of the hermitian models: M.R. Douglas, N. Seiberg and 
S.H. Shenker, Phys. Lett. B244 (1990) 381. 

22. H. Neuberger, Nucl. Phys. B352 (1991) 689. 

23. V. Periwal and D. Shevitz, Phys. Rev. Lett. 64 (1990) 1326; Nucl. Phys. B344 
(1990) 731. 

24. C. Crnkovic, M. Douglas and G. Moore, Nucl. Phys. B360 (1990) 507. 

25. V.G. Kac and M. Wakimoto, Proc. of Symposia in Pure Math., 49 (1989) 191. 



25 



